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Molecular quantum electrodynamics leads to photon frequency shifts and thus to changes in 
condensed matter free energies often called the Casimir effect. Strong quantum electrodynamic 
coupling between radiation and molecular motions can lead to an instability beyond which one 
or more photon oscillators undergo a displacement phase transition. The phase boundary of the 
transition can be located by a Casimir free energy instability. 
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I. INTRODUCTION 



II. PHOTON SCATTERING OPERATOR 



The interaction between electromagnetic field and the 
molecules of condensed matter shifts the frequency of 
the photon oscillators which changes the free energy of 
the condensed matter system. This free energy shift is 
known as the Casimir effect jj]- From a quantum field 
theoretical viewpoint, the Casimir free energy is the con- 
tribution due to one photon loop Feynman diagrams 0- 
Originally, Casimir intended to describe either photon 
oscillator induced forces between a few atoms or between 
fixed perfect conductors. Presently Casimir effects in- 
clude the dielectric and paramagnetic materials as well 
as irreversible effects in moving conductors 0, IE IS- 
The literature is quite extensive but several excellent re- 
views are available 0, IE 0|- Our purpose is to point out 
that if the coupling between the photon oscillators and 
the molecules is sufficiently strong, then the photon os- 
cillators undergo a displacement into a coherent electro- 
magnetic radiation state. While the single photon loops 
do not describe the nature of the radiation coherence, 
the phase boundary between the normal radiation state 
and the coherent radiation state can be located as a free 
energy instability. As the thermodynamic equation of 
state boundary of the normal photon oscillator regime 
is approached, one or more of the photon oscillator fre- 
quencies tends toward zero. Such zero frequency modes 
lead to a Casimir free energy signature of the forthcoming 
radiation transition. 

In Sec^ the connection between the change in the 
photon density of states free energy and the photon scat- 
tering phase shift is reviewed. The Casimir free energy 
then follows by directly superimposing the Planck free 
energies of the shifted oscillators. In Sec lIIII the free en- 
ergy is expressed directly in terms of the (gauge invari- 
ant) electric field photon propagator G and the dielec- 
tric susceptibility x- The Casimir instability is explored 
in Sec lIVI (eventually in terms of the dielectric response 
function e). The phase regime of the Casimir photon os- 
cillator stability is shown to be identical to t he p hase 
regime of the stable Clausius-Mossotti theory [lOl 
The general physical principles involved are discussed in 
the concluding SeclVl 



Consider a condensed matter object which can elasti- 
cally scatter photons. Mathematically the scattering of 
a photon having frequency lj would be described by a 
unitary scattering operator Sioj) whose eigenvalues de- 
termine a set of phase shifts {5y{Lo)} via 



(1) 



The Wigner time delay [13 t,{w) in a photon scattering 
channel v represents the extra amount of time the photon 
remains in the neighborhood of the target due to the 
photon-target interaction. It is well known to obey 



d5^{uj) 
duj 



(2) 



If one sums the time delays over all the channels, then one 
may compute the extra number dN{uj) of photon states 
in a bandwidth diu due the target-photon interaction 



dN(uj) — — Ty{u})du) = — d5,y{ijj). 



(3) 



The Casimir contribution to the condensed matter object 
free energy may be written in terms of the free energy 
f(u!,T) of a simple harmonic oscillator; i.e. 



f{io,T) = -ksTlnl ^ 
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71, n = 0,±1, ±2, Zdois . (4) 



The Casimir free energy is found by simply summing os- 
cillator free energies over all the added photon modes as 
dictated by Eqs.© and It is 



AF = / f{LO,T)dN{u:) 
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(5) 



Now consider the analytic continuation from real fre- 
quency Lo to the complex frequency C in the upper half 
plane 



Q — uj + is where s = '<smC, > 0. 



(6) 



The analytically continued scattering operator obeys the 
dispersion relation 



lndetS'(C^ 



lndet5'(C^) 



w3mlndet5(a;2 + i{)+)dw 

Cj2 - (-2 



(7) 



Our final expression for the Casimir free energy in terms 
of the analytically continued scattering operator follows 
from Eqs.lO) and ((Tjl to be 



AF = ~^ V \ndetS{~u;l). (8) 



4 ^ 
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Let us now consider the more general case wherein the 
scattering operator may also describe inelastic processes 
such as photon absorption by the condensed matter sam- 
ple. The eigenvalue problem, S = e^^^^li/), now 
yields complex phase shifts with 3m > 0. However, 
the Casimir free energy Eq.© remains in tact. This is 
most easily understood in terms of the photon propaga- 
tor. 



III. PHOTON PROPAGATORS 



Starting from Maxwell's equations 
1 fdB\ 

1 fdE\ 4tt^ 



divE = Anp and diwB = 0, 



(9) 



one finds that 



Ml) - A E.-4Jlf + grad.^ (10) 



For any locally conserved charge (dp/dt) + divJ — 0, 
there exists a polarization P such that 



p = —divP and J 



dP 
'dt 



If we invoke the retarded tensor propagator 

2 I ^iQr j c 



G(r,C) = 



for 3m C > 0, (12) 



then Eas. (|10|) and are formally solved by 
E(r,t) = E„(r,t) + 



G r-r',C = * 



P(r',i)dV, (13) 



wherein the incoming photon is described by an electric 
field obeying the vacuum wave equation 



4 (^^ ) - A ^ E„ = 0, 



(14) 



and the outgoing photon has been scattered by the polar- 
ization P. The use of the formal operator replacement 
C — i{d/dt) may be illustrated as follows: The polar- 
ization response in the condensed matter target may be 
related to the electric field by employing a retarded non- 
local response function 

P(r,i) ^ J K(r,r',s) •E(r',i-s)d3r'ds, 



P(r,t) 



K(r,r',s)e-"(^/^*)ds 
d 



E(r',t)d3r', 



P(r,i) = yx(^r,r^C = ^^j •E(r',i)dV, (15) 

with the complex frequency dependent susceptibility de- 
fined as 



x(r,r',C)- / K{r,r', s)e<'ds. 
Jo 



(16) 



The scattering equation follows from Eas. H13() and l|15() 
to be 



E(r,i) = E,„(r,t) 



X r',r",C 



G ( r — r', C = j 



dt 



d_ 

E{r",t)d^r'd^r". (17) 



The Fredholm-Jost determinant [T^ J{() and its log- 
arithmic determinant $(C) for the scattering integral 
Ea. (|17|) are defined as 

J(C) = det{l-G(C)x(C)}, 

$(0 = -lnJ(C) = -rrln{l-G(C)x(C)}- (18) 

Ea. (|18|) is related to the determinant of the scattering 
operator according to 



detS'(C2) = - 



JiC) 



Ji-C) 



(19) 



(11) From Eqs.©, ^ and ^ it then follows that 
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oo ^ 
k=l 



(20) 



Eas. H2Q|l yield the "one loop" contributions to the free 
energy which constitutes the Casimir effect. The Feyn- 
man diagrams are shown in Fig^ If one employs the 
eigenvalue equation for a photon state 



then 



G(C)x(C)k)-^.(C)l^) 



MO" 



(21) 



^(c)-EE^ = -EiMi-^^(c) ■ (22) 



I' k=l 



For the Casimir free energy AF in Eas. lf^ and to 
be stable, it is necessary for the eigenvalues to obey the 
condition 



r7i^(i|ci;„|) < 1 for all v and n, (Stable). 



(23) 



If any eigenvalue (say as a function of temperature) 
passes through unity, then the Casimir free energy will 
have an imaginary part. The transition rate per unit time 
for the system to make a transition into a stable state is 
given by r = — (2/?i)3m AF. The Casimir one photon 
loop free energy in Fig^ is not sufficient for comput- 
ing the free energy of the stable ordered radiative phase. 
Nevertheless, by examining the eigenvalue Eq.j^U equa- 
tion for the mode, i.e. channel which becomes un- 
stable, one may gain insights into the radiative coherent 
state which will be stabilized. The phase diagram for the 
disordered to ordered radiative state may also be com- 
puted. Let us illustrate how this comes about. 



IV. POLARIZATION INSTABILITY 

From Eas. H13|l . (|15|l and H21|l one may examine the 
electric field of a channel mode according to 



G(r-r',C) • x(r',r",C) 

• E,(r")d=^r'd3r" = 77,(C)E,(r)(24) 



Ea. (|24() can be decomposed into two parts: (i) The mode 
E,^(r) produces a polarization Pi/(r). (ii) The polariza- 
tion Pi/(r) radiates the electric field Ej^(r). 

j x(r,r',C)-E,(r')dV = r;,(C)P.(r), 
j G(r-r',C)-P.(r')dV = E,(r). (25) 

If the eigenvalue r] passes through unity, then the field 
mode can be sustained self consistently indicating a ra- 
diative instability. Otherwise, the coupling is two weak 
to maintain the mode and the Casimir free energy con- 
tribution is stable. 

To see what is involved, let us consider a translational 
invariant fiuid for which the bulk Casimir free energy per 
unit volume A/ is determined by electric susceptibility 



x(r-r',C)= / Xk(C) 



?k- (r— r 



(2^ 



(26) 



In detail, Eqs.JEl,®, 123), and (gg) imply 

E(r;k) Eke'''•^ 
Gk(C) - j G(r,C)e-*''"d3r 



-47r 



kk+(C/c)'l 

F - (C/c)2 

?7k(C)Ek = Gk(C)-Xk(C)-Ek, 



2A/ 



= - E yin(l-^k(^K|))^.(27) 



The static (a;„ = 0) contribution to the free energy per 
unit volume in Ea. (|27|l is then determined by 



Afsta 



e(k) = l-r7k(0) 



Ine(k) 



d^k 

(2^ ' 



£(k) = i+An( ^-''if^-^ ) = i+A^mm 



Let us examine Eas. (|12|) and 125|l in the limit of zero 
frequency; i.e. 





FIG. 1: Shown are the one loop Feynman diagrams contribut- 
ing to the Casimir free energy in Eq. 1201 . For each propagator 
wavy line one inserts Gij{ri, rj,Q in accordance with Ea. H12^ . 
For each "self energy part" one inserts Xki{i^k,vi, Q in accor- 
dance with Egs. 11511 and 1161 . The trace "Tr" includes sums 
over polarization indices and integrals over space. 



D(r) 



limG(r-r',C) = VV 
C^o \|r-r' 

linr G(r - r',C) = D(r - r') 



47r 



l^(r-r'), 



3rr — r 1 



Ek(r) = / D(r - r') • Pk{r')d'r' 



(29) 
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A static polarization will produce a local dipolar electric 
field F(r; k) via 

F(r;k) = y D(r-r') •P(r';k)d^r', 

F(r;k) = E(r;k)+ (^^^P(r;k). (30) 

For a fluid with a volume per molecule of v and a nonlocal 
molecular polarizability /3(k), deflned by the local field 



,;P(r;k) =/?(k)F(r;k), (31) 
Eas. H3Q|l and H31(l imply the dielectric response function 



D(r;k) ^ E(r;k) +47rP(r;k) = e(k)E(r;k), 
l + {87r/3(k)/3w} 



e(k) 



1 - {47r/3(k)/3?;} 



(32) 



Note, in the limit of zero wave vector Ea. (|32|l yields 
the usual Clausius-Mossotti formula relating the static 
dielectric constant to the molecular polarizability 

^lim /3(k) — a (polarizability), 

lim £(k) = £ (dielectric constant), 
IkHo 



3v + Sna 



(33) 



3t; — ina 

The Clausius-Mossotti stability condition fTsl IT^ . 

3v > Ana (stable normal phase), (34) 

has then been shown to follow from the Casimir effect 
stability (77 < 1). 

V. CONCLUSIONS 

The quantum electrodynamic interaction between 
molecules and photons gives rise to frequency Lamb 



shiftsflT^ which in turn contribute to a shift in the free 
energy of a condensed matter system. This Casimir con- 
tribution to the free energy can become unstable if the 
interaction coupling is too strong. The instability is the 
signature to a phase transition into a state with coherent 
radiation. We have exhibited such an instability within 
the above statistical thermodynamic considerations. We 
conclude this work by exhibiting the instability directly 
in terms of a collective Lamb shift. 

For zero wave number (k = 0) but finite frequency 
(w ^ 0), the dielectric response function of Eg. 1)32(1 reads 



1 + {87r/3(tj)/3w} 
1 - {ATil3{uj)/iv} 



(35) 



If the molecular polarizability P{lo) can be described by 
a single excitation frequency wq = {AE/h) via 



/3(a;) = 



2 2 ; ' 



then the dielectric response function 



£(C.) 



3u(wo — 0^2) — A-Kuj^a 



(36) 



(37) 



gives rise to a shifted frequency 



1 - — — \ (collective Lamb Shift). (38) 

6V 



The Eq. (|34() for the stable phase regime requires that the 
collective Lamb shifted frequency be real. A strong 
coupling imaginary frequency is a clear indicator of an 
unstable radiation phase. 

While the Casimir effect expressions are sufficient to 
derive the stable photon oscillator regime in temperature 
and density, the one photon loop approximation lacks the 
sensitivity to decide the nature of the radiation ordered 
state. Ordered states of matter are not easily understood 
from a perturbation theory viewpoint. 
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